Landau-Zener transitions in a linear chain. 
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Abstract 

We present an exact asymptotic solution for electron transition amplitudes in an infinite linear 
chain driven by an external time-dependent electric field. This solution extends the Landau-Zener 
theory for the case of an infinite number of states in the discrete spectrum. In addition to the 
transition amplitudes we calculate the effective diffusion constant. 
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1 



Landau-Zener (LZ) theory!,! treats a quantum system placed in a slowly varying external 
field. If such a system was prepared in a state of its discrete spectrum, it adiabatically follows 
this state until its time dependent energy level crosses another one. Near the crossing point 
the adiabaticity can be violated and the system can escape from the state it occupied initially 
to another one. Landau and Zener found the transition probability for two-level crossing. 
The crossing of more than two levels at the same time is generally an unlikely coincidence. 
However, in some systems such a multi-level crossing may occur systematically, due to the 
high symmetry of the underlying Hamiltonian. The transition matrix for special cases of 
multi-level crossing was studied in Refs.S&iffl!. Presently only a few exact results for multi- 
level crossing are known. One of them relates to a multiplet of atomic electronic states with 
a total spin S or total rotational moment J larger than 1/2 in a varying external magnetic 
fieldH. The Zeeman splitting between 25* + 1 or 2 J + 1 levels regularly vanishes at nodes 
of the magnetic field. Another exactly solvable model displaying multi-level crossing is the 
so-called bow-tie model!, whose physical interpretation is not obvious. 

Since its creation in 1932, LZ theory has had numerous applications. They include 
molecular pre-dissociation0, slow atomic and molecular collisions^, and electron transfer 
in biomoleculesS Recently Wernsdorfer et aZ.0,0 employed the LZ theory to describe 
consistently the step-like shape of the hysteresis loop in special molecules with large magnetic 
moments called nano-magnets. Using the LZ probability formula these authors were able to 
find the extremely small tunnel splitting of the classic degenerate ground states and even 
to reveal oscillations of this value in external magnetic field. This beautiful experiment, 
together with its clever treatment is a new triumph of quantum mechanics and, in particular 
LZ theory. 

The problem considered in this article is closely related to another application of LZ 
theory: electronic transfer in donor-acceptor complexes!. In this process of biological and 
chemical importance, an electron tunnels between initial and final positions through a long 
chain of identical sites. There are two limiting cases for such a process. In the first case 
there is no coherence between two sequential tunnelling processes connecting nearest neigh- 
bor sites. In this case the probability of tunnelling through several sites is very small in 
comparison to that for one-site tunnelling. This limiting case was studied earlier!. We 
consider the opposite limiting case in which the sequential tunnelling processes are highly 
coherent and tunnelling through many sites becomes available. 
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If the coherence between LZ transitions is lost, the problem is reduced to multiplication 
and addition of probabilities, each described by a proper LZ expression. The price we must 
pay for incorporating the coherence between different transitions is a strong reduction of 
the class of quantum systems considered. The number of crossing levels in such systems 
must be infinite. The hopping amplitudes from a site to its neighbors must be all identical. 
Physically it describes the quantum electron transfer between donor and acceptor separated 
by a long polymer strand (molecular bridge). The bridge can be considered as a linear 
array of identical sites. Such one-dimensional atomic-scale wires were intensely studied, 
both experimentally and theoreticallyEl'Ei0. Our result can be also applied to transitions 
among electron states in semiconductor superlatticesSl!. 

We study the tunnelling of a particle in such systems driven by a time- dependent ho- 
mogeneous external field. An important assumption is that all molecular fragments in the 
chain are identical. An electric field splits the energy levels at different sites of the chain 
and suppresses the transitions, which occur within a narrow intervals about times when the 
electric field becomes zero. Since the tunnelling is a fast process, we disregard the oscillatory 
relaxation originating from phonons and other elementary excitations. 

Let denote |n > a state located at the n-th site of the chain. We assume that these 
states form a complete orthonormal set (Wannier basis). In terms of this set the electron 
Hamiltonian reads: 



N 

H = I n)(n + l | +c.c.)+F(t)n | n)(n |; (1) 

n=l 

F(t) = eE(t)a 

where E(t) is the electric field, e is the electron charge, a is the distance between sites and 
7 is the coupling constant (hopping amplitude). A series of exact solutions for the time- 
dependent Shrodinger equation with the Hamiltonian (]]]) for N = oo, known as drifting 
plane waves, was found long Below we solve the same problem for fixed 

initial conditions, thus resolving the multi-state LZ problem. 

The states |n > are conventionally called the diabatic states. They are the eigenstates 
of the diagonal part of this Hamiltonian Hq = F(t)n | n)(n |. The eigenstates of the 
total Hamiltonian ([!]) depending on t or F(t) as parameters are called adiabatic states. 
Until | F(t) 7, the diabatic levels are close to the adiabatic ones and the transitions 
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FIG. 1: Dependence of eigenvalues of the Hamiltonian (Q) on E. Solid lines show adiabatic energies 
as function of E, dashed lines depict diabatic energies. While the diabatic levels intersect in one 
point, the adiabatic levels do not intersect. 

between levels are suppressed. This is the adiabatic regime. The adiabaticity is violated 
in the vicinity of the electric field nodes determined by the inequality | F(t) |< 7, where 
all transitions proceed. By level crossing we mean that the diabatic levels cross, the exact 
eigenvalues of the Hamiltonian ([!]) never cross, in accordance to the Wigner - Neumann 
theorem. It is convenient to place the time origin t — directly at the node of F(T). Since 
only a narrow vicinity of the node is substantial for transitions, the exact dependence of the 
field on time can be reasonably approximated by linear one: F(t) ~ -^(0) -t. At zero electric 
field E and free boundary conditions the Hamiltonian ([!]) can be diagonalized analytically. 
Its spectrum is: 

e j = 2gcoa(Trj/N); j = l...N (2) 

For non-zero field we have found the adiabatic eigenvalues numerically. The result is 
shown in Fig. 1 for a finite chain with 15 sites. For comparison the diabatic levels are 
depicted in the same figure. 

We proceed to solve the time-dependent Schrodinger equation with the Hamiltonian ([[]) . 
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Its matrix representation reads: 



H nm = nF(0)t5 nm + 7(5 m , n +i (3) 

For an infinite chain (N — > oo) equation (|3|) is valid for all n and m. After a proper rescaling 
of time the Hamiltonian (Q) becomes dimensionless: 

(4) 



It depends on only one dimensionless number, g = j/yF(0), which is the Landau-Zener 
parameter. Let the time-dependent state vector be | a,t) — I] n c n (t) | n). Then the system 
of equations for the amplitudes c n (t) reads: 

ic n = ntc n + g(c n -i + c n+ i) (5) 

The transition matrix element T Uin i should be identified with the t — ■> +oo asymptote of an 
amplitude c n (t) for a solution obeying the initial condition |c m (t)| 2 = 5 m y at t — > — oo. Since 
all c m (t) except of c' n (t) are zero at t — > — oo, the initial condition can be more explicitly 
written as 

c m (t -> -oo) = 8 m>n/ exp (-in't 2 /2) (6) 

We multiply the asymptotic values of c n (t) by exp [int 2 /2] to remove strongly oscillating 
phase factors from T n y. 

Now introduce an auxiliary function u(ip,t) = J2^=-oo c nif)^ mip ■ The system is equiv- 
alent to the following equation in partial derivatives for u(ip,t): 

du du 

~dt +t d^ + 2l9UCOSV = ^ 

The initial condition (|6|) is equivalent to the initial condition: u(<p, t — > — oo) — > exp[m'(— ^ + 
ip)]. Given the solution u(ip,t), the amplitudes c n (t) can be found by the inverse Fourier 
transformation c n (t) = j- J 27r u((p,t)e~ mip d(p. The solution of eq. (0) that obeys proper 
boundary conditions is: 

/ * 

u((p, t) = exp 



-i[2g I cos Up - l — + y j dt' + 



'( ^ 



(8) 



Putting t = +00 in the solution @ and taking the inverse Fourier-transform, we arrive at 
following asymptotic values: 

/ int 2 (n' — n)7r\ ,„ r— . 

c n (t) « exp I Y + 4~J J \n-n'\(zV^g) (9) 

Thus, the scattering amplitudes in terms of modified states, with the fast phase factor 
exp(— int 2 /2) incorporated, are: 



T n , n , = (n\T\n')= e^'-^J^^V^g), (10) 

where the operator T is expressed in terms of the evolution operator U (t, if) for the Hamil- 
tonian (|J) in the interaction representation: 

T= lim exp(i [ H (T)dT)U(t,t')exp(i F H (t)(1t) (11) 

t^oo;t'^-oo Jt> JO 

The matrix elements T n>n > display an infinite number of oscillations with the LZ parameter 
g. However, for large \n — n'\ the oscillations start with g > \n — n'\. These oscillations 
can be observed experimentally by varying the field sweep rate E(0). For small values of g 
the amplitudes are small and quickly decrease with growing \n — n'\. In Fig. 2 we depict 
transition probabilities for several levels n closest to the initial one nf versus the Landau- 
Zener parameter g. Fig. 3 shows the dependence of the transition amplitude on | n — n'\ at 
a fixed value of g. 

For large g ^> \n — n'\ the asymptotic values of the amplitudes (|T0|) are: 

{n 1 1 - > ~ jym^ cos(2 ^ s " ~ ?' (12) 

It is instructive to compare this result with other exactly solvable generalized Landau- 
Zener models. Most of them refer to systems with a finite number of states N. In the 
limit g ^> N the transition probabilities behave like an exponent exp(— C(n, n')g 2 ), where 
the C(n,n') do not depend on g. In contrast, the result ( J12] ) displays a power law with 
oscillations instead of an exponential dependence on g for large g. This is the manifestation 
of quantum interference of different Feynman trajectories, which are discrete in the chain. 
A step in a trajectory has average length g (see below). Such a step cannot be realized in a 
system with a finite number of states if g ^> N. 

The mean square displacement at one crossing event is: 
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FIG. 3: Transition probabilities vs. n — n' at a fixed g = 1.6. Solidline is a guide for eye. 
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< (n-n') 2 >= £ (n-n') 2 \J ln ^ nll (2V2^g)\ 2 = 

n=— oo 

= 47H7 2 (13) 

If the external field is periodic in time and the coherence between crossing events is lost!, the 
electron performs a random walk, i.e. it diffuses. Assume the field to oscillate harmonically 
as F(t) = F sm(ut). At the nodes tj, = irk/u (k is an integer) all diabatic levels cross 
together. The squared Landau-Zener parameter is g 2 = -j*— . The diffusion coefficient is 



,2 ^ f„ m t\2 



D = 2a < {n — n') > /T, where T = 2tt/lu is the period of oscillations and the factor 2 
accounts for two crossing events per period. Collecting these results and equation (|T3D , we 
find: 

D - ^ (14) 

This result does not depend on the frequency of the external field. 

The theory can be extended to a more general Hamiltonian incorporating hopping be- 
tween any two sites, but conserving translational invariance: 



H = ^2 H rn , n \m) (n\; 

m,n 

H mn = TltS mn + g\m- n \] 9—k = 9k (15) 

For simplicity we present below the result for real hopping amplitudes gk = g~k- 



J(n'-n)Tr/4 2n 



( n |T n J = — J exp (-i2V2irf((p,gj) 

+i(n — njip) dip (16) 



o 



where f(tp, gj) =Ek^ cos kip. 

The model (|1|) can be generalized also to incorporate internal degrees of freedom of 
identical chain fragments. In this case the local states are described by amplitudes a nj0l with 
two indices. The first index n denotes the position and the second index a labels the inner 
states. The Schrodinger equation for the amplitudes then reads: 



ia n , a = (nt + e a + 5 a t)a n , a + ^ g a A a n+i,p + a>n-i,p) (17) 



S 



where the indexes a, (3 = 1 . . . Ni nt run over the internal states of the molecular wire segment. 
Changing to variables a nj0l = b n)0t e~ int2 1 2 we eliminate the term proportional to t in equation 
(pp. Introducing a new function u a (ip) = J2^=-oo^n,a^ nip , we reduce the infinite system 
(PI) to a finite set of N int ordinary differential equations: 

iu a = (e a + 5 a t)u a + 2 cos(t 2 + (p) ^ g a ,pup (18) 

P 

in which <p plays the role of a parameter. The initial conditions are b n ^ a (t — > — oo) = 
exp (— id a t 2 /2)S n y and u a ((p,t — > — oo) = exp [(— iS a t 2 /2) +in'(p]. Thus, the variable (pa- 
rameter) ip enters not only in the system flUf), but also in the initial conditions. This 
system must be solved for all values of parameter <p in the interval (0, 2tt) . The inverse 
Fourier-transformation yields the evolution operator just as for the case N int = 1. An 
analytical solution of the system fll8|) is possible for some special choices of the parame- 
ters e a , 5 a and g ay p. For example, two identical coupled chains correspond to Ni nt = 2. 
Then the indices a, (3 take on the values 1,2. The simplest solvable choice of parameters is: 
E\ = 0,e 2 = £; #i = 0, 5 2 = S',gi,i = 92,2 = 7! g\,i = ^2,1 = l'- Exact solution of this model 
can be reduced to ones solved in this article together with solvable two level LZ model. 

In conclusion, we have generalized the LZ theory to an infinite number of crossing levels. 
Physically it describes an electron on an infinite chain subject to a time- dependent electric 
field. The high symmetry of the problem allows us to find not only the asymptotics, but also 
the intermediate values of the amplitudes. Our solution is valid for an infinite chain with 
translational symmetry. We demonstrated it for the case of a simple primitive cell, but it 
can be generalized and in some cases exactly solved even if the primitive cell contains more 
than one site. 

Finite size effects do not permit us to apply our result ( |10|) directly to the transition 
amplitude from the first to the last site of the chain even if it contains many sites. In 
particular, the transition amplitude from one end to another does not oscillate. However, 
our calculation of the diffusion coefficient (0) is valid since diffusion presumably proceeds 
far from the ends of the chain. Certainly, we assume that coherence is lost during the time 
interval between two sequential crossing events. 

Our solution demonstrates a phenomenon that is probably common for most systems 
with multi-level crossing: oscillations of the transition probabilities as a function of the LZ 
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parameter and site position (distance between diabatic levels). However, their asymptotic 
values for large values of the LZ parameter differ from those for other solvable multi-state 
LZ-models with a finite number of states. We expect that in a general situation with N 3> 1 
crossing levels, the transition probabilities will behave similarly to those found in this work 
for 1 ^ g ^ N, provided that the initially occupied states are far enough from the diabatic 
spectrum boundaries. 

Finally we discuss the relationship between our problem and a typical problem for semi- 
conductor superlattices@'Ei. The latter is associated with Anderson localization. The di- 
abatic levels at sites are randomly distributed. In one and two dimensions all sites are 
localized. If the width of the energy distribution A is much less than the tunnelling am- 
plitude 7, the localization length in one dimension is 07/ A. To enhance the tunnelling 
through a chain it is reasonable to apply a time-dependent electric field. The electric field 
is substantial if F7/A > A where F is a typical value of F(t). Our approximation is valid 
if the inequality is strong: F7/A ^> A. Tunnelling transitions in the field proceed during 
an interval of time r defined by relation F(r) ~ 7. The value F(r) can be accepted for F . 
We see that the strong inequality 7 ^> A guarantees the existence of the strong field limit 
in which the randomness of levels can be ignored. This requirement does not impose any 
limitations on the LZ parameter g. 
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